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Results are reported for N* masses in the 1/2+ and 1/2— sectors on an anisotropic lattice. The gauge action 
is the usual plaquette plus rectangle type, the quark action is of the D234 type, both having tadpole-improved 
tree-level coefficients. Clear splittings from the nucleon ground state are observed with smeared operators and 
500 configurations. The first excited states in each channel are further isolated, using a correlation matrix and 
the variational method. The basic pattern of these splittings is consistent with experimental observations. 



1. Introduction 

There is increasing experimental information 
on the N* spectrum from JLab and other acceler- 
ators, and the associated desire to understand it 
from first principles. Given that state-of-the-art 
lattice QCD simulations have produced a ground- 
state spectrum that is very close to the observed 
values jjj, it is important to extend beyond the 
ground states. The rich structure of the N* spec- 
trum j| provides a fertile ground for exploring 
how the internal degrees of freedom in the nu- 
cleon are excited and how QCD works in a wider 
context. One outstanding example is the splitting 
pattern between the observed 1/2+ spectrum and 
its parity partner: the 1/2— spectrum. The split- 
tings are a direct manifestation of spontaneous 
chiral symmetry breaking of QCD, because with- 
out it QCD predicts parity doubling in the baryon 
spectrum. 

Several lattice studies of the N* spectrum 
have appeared recently. In H the 1/2+ spec- 
trum was first explored with a 0(a 2 ) improved 
next-nearest-neighbor action. In [Q the splitting 
N1/2- — N1/2+ was examined with the domain- 
wall fermion action. In [gj results are reported 
for the Ni/ 2 ± and the A 3 / 2 ± based on the non- 
perturbative 0(a) improved clover action. All the 
calculations show a clear splitting of the N 1 / 2 - 
from the ground-state nucleon. In this work, we 



report further calculations in the 1/2+ sectors us- 
ing a highly-improved action on a anisotropic lat- 
tice. 

2. Method 

Two independent, local interpolating fields 
coupling to the nucleon are considered: 

Xi(x) = e abc (u Ta (x)C l5 d b (x)) u%x) , (1) 

X2 ( x ) = e abc (u Ta (x)Cd b (x)) l5 u c (x) . (2) 

The interpolating fields for other members of the 
octet can be found by appropriate substitutions of 
quark fields |J. Despite having explicit positive- 
parity by construction, these interpolating fields 
couple to both positive and negative parity states. 
A parity projection is needed to separate the two. 
In the large Euclidean time limit, the correlator 
with Dirichlct boundary condition in the time di- 
rection and zero spatial momentum becomes 



G(t) = ^2<0\ X (x)x(0)\0> 



X 



2 74 + 1 M + t 

2 



+ A 2 



-74 



-M- t 



(3) 



(4) 



"Lattice Hadron Physics Collaboration 



The relative sign in front of 74 provides the so- 
lution: taking the trace of G(t) with (1 ± 74)/4 
respectively. A method on how to do the projec- 
tion at finite momentum was also proposed in j|] 
which involves only diagonal elements. Three 
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types of correlators are considered in this work: 
two diagonal ones Ni =< X1X1 > and N% =< 
X2X2 >, and the crossed iV 3 =< X2X1 > + < 
X1X2 > 

The anisotropic gauge action of m, and the 
anisotropic D234 quark action of || are used. 
Both have tadpole-improved tree-level coeffi- 
cients. Similar anisotropic actions have been 
used to study systems with heavy quarks jo|— pTT| . 
One advantage of anisotropic lattices is that with 
modest lattice sizes one can access large spa- 
tial volumes while having a fine temporal resolu- 
tion, important when the states are extended and 
heavy. A 10 3 x 30 lattice with a s « 0.24 fm and 
anisotropy £ = a s /at — 3 is used. The tadpole 
factor from the plaquette mean link is u s — 0.81 
while ut is set to 1. In all, 500 configurations are 
analyzed. On each configuration 9 quark prop- 
agators are computed using a multi-mass solver, 
with quark masses ranging from approximately 
780 to 90 MeV. The corresponding mass ratio 
ir/p is from 0.93 to 0.60. The source is located 
at the 2nd time slice. A gaussian-shaped, gauge- 
invariant smearing function fl2|| was applied at 
both the source and the sink to increase the over- 
lap with the states in question. Statistical errors 
are derived from a bootstrap procedure. 

3. Results and Discussion 

Figure [l] shows some sample effective mass 
plots at quark mass no. 3 and 7 for all three corre- 
lators in the source-sink combination of SL which 
has the best signal among LL,LS,SL,SS. Valid 
plateaus exist in all cases. Figure |2| shows the ex- 
tracted masses as a function of the quark mass. 
The Ni and N3 give roughly the same mass, while 
N 2 yields a consistently higher mass. It turns out 
that xi has very little overlap with the ground- 
state nucleon and couples mostly to the first ex- 
cited state, the Roper A r 1 / 2 +(1440). This point 
will be clarified below. All three negative-parity 
correlators give consistent masses. The ordering 
of A r i /2+ (1440) and AT 1/2 _(1535) at the quark 
masses considered is inverted compared to that 
from experiment. These results are consistent 
with those observed in Q. It would be interest- 
ing to see if the ordering is reversed at smaller 
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Figure 1. Effective mass plots at two quark 
masses. 
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Figure 2. Masses of N x / 2 ± as a function of quark 
mass for all three types of correlators. 
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quark masses where meson cloud becomes more 
important. The crossing can be expected even in 
the quenched approximation where quarks wig- 
gling backward in time provide part of the meson 
cloud. Of course quenched QCD is sick at small 
quark masses, so a correct description must in- 
clude the sea quarks. There is already evidence 
of such crossing at very small quark masses in the 
case of ao and a\ mesons H@. The effects of 
the meson cloud show up as non-analytic behav- 
ior in the quark mass. Some curvature is already 
apparent in the data so a linear extrapolation is 
undesirable. 

In order to make contact with experiment, Fig- 
ure H presents mass ratios of the 1/2— states to 
the nucleon ground state as a function of (w/p) 2 . 
This way of presenting data has minimal depen- 
dence on the uncertainties in determining the 
scale and the quark masses. These ratios appear 
headed in the right direction compared to experi- 
ment. The A]/ 2 - (1405) deserves special attention 
since it requires a larger curvature. At present it 
is almost degenerate with Af]/ 2 _(1535) while in 
the physical spectrum its mass is lower despite 
having a heavier strange quark. This has been 
a long-standing puzzle. The origin of the mass 
suppression of the Ai/ 2 _(1405) lies in strong res- 
onance channels such as the KN and 7rS which 
would present themselves as nonlinearities in the 
quark-mass extrapolation. Whether it happens 
at smaller quark masses remains to be seen. 

To extract information on the first excited 
states in each sector, the 2 by 2 correlator matrix 
from the two interpolating fields is considered, 
using the variational method originally proposed 
in |l5j . Figure ^ shows the results in the nucleon 
channel. It confirms that \2 couples mostly to 
the the iV 1/2+ (1440). The splitting of the Roper 
N\/2+ (1440) from the nucleon is roughly 5 times 
that of the N 1/2 _ (1620) from the iV 1/2 _(1535), 
consistent with experimental observations. This 
may be traced back to the structure of the in- 
terpolating fields used to excite the states. The 
Xi and xi have upper spinor components of the 
order 1 and (p/E) 2 respectively, whereas the in- 
terpolating fields for 1/2- states, Xi = 75Xi an d 
X2 = 75X2, both have upper spinor components 
of the order p/E which also introduce deriva- 
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Figure 3. Mass ratio N*/N as a function of 
(ir/p) 2 . The experiment points are indicated on 
the left. 
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Figure 4. Ground states and 1st excited states in 
the nucleon channel from the correlator matrix. 
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tives of the quark field operators. In the quark 
model, the two nearby low-lying N 1 / 2 - states of 
the physical spectrum are described by coupling 
three quarks to a spin-1/2 or spin-3/2 spin-flavor 
wave function coupled in turn to one unit of or- 
bital angular momentum. In a relativistic field 
theory one does not expect Xi an d X2 to isolate 
individual states. However, it is expected that Xi 
will predominantly excite the lower-lying j = 1/2 
state associated with I = 1 and s = 1/2 while X2 
will predominantly excite the higher-lying state 
associated with s = 3/2. This expectation is 
borne out in the variational analysis. The results 
should improve with a larger basis of interpolat- 
ing fields. An independent method based on max- 
imum entropy |l6| ] is being applied to verify the 
results. 

4. Conclusion 

The potential of anisotropic lattices, smeared 
operators and correlator matrix is demonstrated 
in probing the low-lying N*s in the 1/2+ and 
1/2— sectors. The pattern of the splittings is 
mostly consistent with experiment. More defi- 
nite comparisons should address the systematics: 
finite a, finite volume, statistics, and most im- 
portantly, the chiral limit. There are still unre- 
solved issues, such as the inverted ordering of the 
iV 1/2+ (1440) and the A 1/2 _(1535), and the de- 
generacy of the Ai/2-(1405) and the ATi/2_(1535) 
at present quark masses, that depend on this 
limit. One way to proceed is to use ChPT to ex- 
trapolate the lattice data at quark masses where 
the quenched approximation is benign, like that 
proposed in JljJ . The more desirable way is to be 
able to simulate directly at or near the physical 
quark masses. Toward this end, good chiral prop- 
erties are essential, and the domain-wall and the 
overlap quarks [|l8|,[l9| seem to hold the promise. 
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